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ABSTRACT 

Observations of the universe reveal that matter within it clusters on a variety of scales. 
On scales between 10 - 100 Mpc, the universe is spanned by a percolating network of 
superclusters interspersed with large and almost empty regions - voids. This paper, 
the first in a series, presents a new ansatz which can successfully be used to determine 
the morphological properties of the supercluster-void network. The ansatz is based on 
a surface modelling scheme (SURFGEN), developed explicitly for the purpose, which 
generates a triangulated surface from a discrete data set representing (say) the dis- 
tribution of galaxies in real (or redshift) space. The triangulated surface describes, at 
progressively lower density thresholds, clusters of galaxies, superclusters of galaxies 
and voids. Four Minkowski functionals (MFs) - surface area, volume, extrinsic curva- 
ture and genus - describe the geometry and topology of the supercluster-void network. 
On a discretised and closed triangulated surface the MFs are determined using SUR- 
FGEN. Ratio's of the Minkowski functionals provide us with an excellent diagnostic 
of three dimensional shapes of clusters, superclusters and voids. Minkowski function- 
als can be studied at different levels of the density contrast and therefore probe the 
morphology of large scale structure on a variety of length scales. Our method for 
determining the Minkowski functionals of a triangulated iso-density surface is tested 
against both simply and multiply connected eikonal surfaces such as triaxial ellipsoids 
and tori. The performance of our code is thereby evaluated using density distribu- 
tions which are pancake-like, filamentary, ribbon-like and spherical. Remarkably, the 
first three Minkowski functionals are computed to better than 1% accuracy while the 
fourth (genus) is known exactly. SURFGEN also gives very accurate results when ap- 
pfied to Gaussian random fields. We apply SURFGEN to study morphology in three 
cosmological models, ACDM, tCDM and SCDM, at the present epoch. Geometrical 
properties of the supercluster-void network are found to be sensitive to the underlying 
cosmological parameter set. For instance, the percolating supercluster in ACDM turns 
out to be more filamentary but topologically simpler than superclusters in rCDM and 
SCDM. It occupies just 0.6% of the total simulation-box volume yet contains about 4% 
of the total mass. Our results indicate that the surface modelling scheme to calculate 
Minkowski functionals is accurate and robust and can successfully be used to quantify 
the topology and morphology of the supercluster-void network in the universe. 



Key words: cosmology: theory — structure formation — morphology — superclusters — 
voids-cosmic web 
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1 INTRODUCTION 

The existence of the supercluster-void network is one of the 
most intriguing observational features of our universe. Red- 
shift surveys of galaxies confirm that on very large scales 
the bulk of matter in the universe is concentrated in clus- 
ters and superclusters of galaxies which are separated by 
large almost empty regions, appropriately called voids. At 
moderate density thresholds (5 ~ 1, the supercluster network 
percolates even though it occupies a tiny fraction of the total 
volume, and so has a small filling fraction. 

The morphology of the supercluster-void network is 
quite complex and has inspired evocative descriptions such 
as being 'honey-comb-like', 'bubble- like', 'a filamentary net- 
work', 'swiss cheese', 'cosmic web' etc. Indeed the intricate 
weave of large scale structure is becoming increasingly evi- 
dent as results from progressively larger galaxy redshift sur- 
veys (CfA, IRAS, PSCz, LCRS, 2dFGRS, SDSS) demon- 
strate. Quantifying the properties of the supercluster-void 
network is clearly one of the cardinal tasks facing cosmol- 
ogy today. Since the network evolved from an almost fea- 
tureless Gaussian random field (as suggested by observa- 
tions of the cosmic microwave background) it is impor- 
tant to understand how such a rich and complex cosmic 
tapestry - characterized by prominent non-Gaussian fea- 
tures (clusters, superclusters, voids) - could have arisen 
from small and statistically random initial conditions. Fur- 
thermore, since the network is highly evolved it is un- 
likely that its principle features can be described within a 
strictly perturbative framework which breaks down when 
5 >, 1. It is therefore both interesting and revealing that 
approximations which successfully describe non-linear dy- 
namics (such as the Zeldovich approximation and the Ad- 
hesion model) generically predict the formation of filaments 
and pancakes, which interweave an d percolate to form the 
large s cale structure of the univers e fShandarin fc Z eldovichI 
Il989l: iGurbatov. Saichev fc Shandarin 1985, 1989^) . Fur- 
thermore, this interweaving network of filaments, sheets 
and voids - predicted both by the Zeldovich approxi- 
mation and the adhesion model - is readily seen in N- 
body simulations of gravitation al clustering for a wide va- 
riety of initial power spectra (iKlxBin, fc ShandarinI Il993t 



ISathvaprakash. Sahni fc Shandarinl ll99d) . It therefore ap- 
pears that a filamentary distribution of large scale structure 
is an almost generic outcome of pressureless (dark) mat- 
ter clustering from Gaussian initial conditions. (A theoret- 
ical model for understanding the morphology of large scale 
structure can be found in the "co smic web" hypothesis of 
iBond. Kofman fc Poeosvanl il99d) . Insightful overviews of 
the statistical and dynamical techni ques used in studies o f 
large scale structu re can be found in:| Sahni fc Cole3jl99El) : 
IMartinez fc Saail ll200lD : lvan de Weveaertl J2002D .') 

A number of statistical measures have been ad- 
vanced to quantify the pattern made by galaxies as 
they cluster in our uni verse. Prominent among them 
are P e rcolation Analysis iZeld ovich, Einasto fc ShandarinI 
Il982l: ISh andarin fc Zeldovich' '1983^. Coun ts in CeUs 
jj^^s"fc D emarg ue 1983; de LaoDaront, GcUe r fc Huchraj 



), Minimal Spanning T rees IB arrow. Sonoda fc Bhavsar 
198E ) , the Genus measure (iGott. Melott fc Dickinsonnl986ll 



etc. (To some extent these methods complement the tra- 
ditional approach to quantify clustering using the hier- 



archy of correlation functions jPeeblej Il980l) . which be- 
come cumbersome to evaluate beyond the three point 
function.) Recently Mecke, Buchert and Wagner (1994) 
have introduced the Minkowski functionals (hereafter 
MFs) to cosmology. MFs contain information pertain- 
ing to geometry, connectivity (percolation) and topol- 
ogy (genus). In addition, the ratio's of Minkowski func- 
tionals quantify the morphology of large scale structure. 
This has led to the construction of a set of measures 
called Shapefinders which tell us whether the distribu- 
tion of matter i n superc lustcrs/voids is spherical, planar, 
filamentary etc. dSahni. Sathvaprakash fc Shandarin 199^; 
ISathvaprakash. Sahni fc Shandarinlll998f) . (Applications of 
Shapefinders to gala xy catalogues and N-body simulations 
has be e n discussed in[B asilakos^ Plionis fc Row an-RobinsonI 
j2000l) : lK'olokot ronis. Basilakos fc Plionis? feOOlV Moment- 
based studies of morphology predating the Shap efinders 
can be found in fBabul fc StarkmanI lll99d) : lLuo fc Vishniad 
il995ll .) Thus between them, the four Minkowski functionals 
contain valuable information regarding both the geometrical 
as well as topological distribution of matter in the universe. 
Since Minkowski functionals are additive in nature one can 
glean information regarding both individual objects (galax- 
ies, clusters, voids) as well as describe the supercluster-void 
network in its totality. 

There have been three major attempts made to study 
the morphology of the LSS using MFs; these efforts differ in 
their approach of evaluating the MFs: (1) Firstly, Boolean 
grain models study the MFs of surfaces which result due to 
intersecting spheres decorating the input point-set (Mecke, 
Buchert fc Wagner 1994). (2) Secondly, Krofton's formulae 
make it possible to calculate MFs on a density field defined 
on a grid (Schmalzing fc Buchert 1997). In this case the 
MFs are calculated by using the information of the number 
of cells in ID (vertices), 2D (faces) and 3D (cuboids). (3) 
Finally, an alternative, resolution-dependent approach con- 
sists in employing the Koenderink invariants (Schmalzing fc 
Buchert 1997, Schmalzing et al. 1999). 

This paper introduces a radically new and power- 
ful approach for computing the Minkowski functionals. 
This approach consists of constructing iso-density sur- 
faces using an elaborate surface modelling scheme de- 
fined in terms of excursion sets of a density field. 
These surfaces are triangulated and the MFs are eval- 
uated for the resulting closed polyhedral surface. A 
similar algorithm has b een used in studies of two- 
dimensional CMB maps fN ovikov. Feldman fc ShandarinI 
^99; Shandarin ct al. 2002). However, developing a three- 
dimensional algorithm requires a new component that builds 
a surface of approximately constant (to linear order) density. 
It has been implemented and reported in this paper. 

The MFs so constructed are tested against known re- 
sults for simply connected surfaces (triaxial ellipsoids) as 
well as multiply connected surfaces (triaxial tori). In all 
cases we find that our ansatz for determining MFs on poly- 
hedra gives results which are in excellent agreement with ex- 
act continuum values for these quantities. We have further 
tested the performance of SURFGEN against Gaussian ran- 
dom fields for which the analytic prediction for MFs is avail- 
able and find that our calculations reproduce the analytic 
results to a remarkable precision. This gives us confidence 
that our ansatz is well suited for application to large N-body 
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simulations and three-dimensional galaxy redshift surveys. 
Encouraged by the success on these preliminary tests, we ap- 
ply our method to simulations of three cosmological models 
due to the Virgo consortium. 

The rest of this paper is arranged as follows: In Sec- 
tion 2, we provide a brief overview of Minkowski function- 
als and present an algorithm which computes MFs for a 
triangulated polyhedral surface. In Section 3, we describe 
the SURFGEN algorithm which triangulates a general class 
of surfaces including iso-density surfaces describing super- 
clusters or voids. In Section 4 we test our method against 
standard eikonal surfaces — spheres, ellipsoids, tori — and 
demonstrate the accuracy of our algorithm. Special bound- 
ary conditions are needed to reproduce the analytic expres- 
sions for MFs of a Gaussian random field. The related dis- 
cussion along with the results for Gaussian random fields is 
presented at the end of the paper in an appendix. Section 5 
is devoted to a detailed morphological study of cosmological 
N-body simulations. The models investigated are ACDM, 
tCDM and SCDM at the present cosmological epoch. Our 
conclusions are presented in Section 6 which also discusses 
possible further applications of SURFGEN and the MFs. 



2 MINKOWSKI FUNCTIONALS AND 
TRIANGULATED SURFACES 

The large scale structure of the universe can be studied 
on various scales by considering the geometry and topol- 
ogy of excursion sets of the density contrast (5(x) defined 
as = {x|5(x) ^ 5th}, for overdense regions (clus- 

ters, superclusters) and = {x|(5(x) ^ Sth}, for under- 
dense regions (voids). By specifying a given density thresh- 
old one effectively defines an iso-density surface over which 
the Minkowski functionals should be evaluated. Depending 
upon whether the surface encloses overdense region of the 
space or underdense region, the surface refers to a clus- 
ter/supercluster or a void. Note that the term 'cluster' is 
used for any region in the excursion set {5 > 5th) which 
in general is different from the astronomical definition. The 
terminology refiects the cluster analysis jargon. The follow- 
ing four Minkowski functionals describe the morphological 
properties of an iso-density surface in three dimensions ^ 

• [l]Area 5* of the surface, 

• [2] Volume V enclosed by the surface, 

• [3]Integrated mean curvature C of the surface (or inte- 
grated extrinsic curvature), 

where Ri and R2 are the principal radii of curvature at a 
given point on the surface. 

• [4] Integrated intrinsic (or Gaussian) curvature x of the 
surface - also called the Euler characteristic 

^ There are (n + 1) Minkowski functionals in n dimensions each 
defined for an (n — 1) dimensional hypersurface. Two dimensional 
MFs have been used to analyze the LCRS slices by Bharadwaj 
et al. (2000) and the anisotropy of the cosmic microwave back- 
ground by Schmalzing & Gorski (1998) and by Novikov, Feldman 
& Shandarin (1999). 




A related measure of topology is the genus G — 1 — x/2- 
^ Multiply connected surfaces have G > while G — for 
a simply connected surface such as a sphere. (The topolog- 
ical properties of all orientable surfaces are equivalent to 
those of a sphere with one or more 'handles'. Thus a torus is 
homeomorphic to a sphere with one handle, while a pretzel 
is homeomorphic to a sphere with two handles etc.) While 
the genus provides information about the connectivity of a 
surface, the remaining three MFs are sensitive to local sur- 
face deformations and hence characterize the geometry and 
shape of large scale structure at varying thresho lds of the 
density dSahni. Sathvaprakash fc ShandarinlllQQST) . 

In nature one seldom comes across surfaces that are 
perfectly smooth and difi^erentiable. As a result expressions 
[1] - [4] which make perfect sense for manifolds C, n ^ 2 
are woefully inadequate when it comes to determining the 
Minkowski functionals for real data sets which are grainy 
and quite often sparse. Below we describe how one can de- 
termine the Minkowski functionals (and derived quantities) 
for surfaces derived from real data by interpolation. The in- 
terpolation scheme which describes iso-density contours as 
triangulated surfaces will be presented in detail in the next 
Section. 

We construct a polyhedral surface using an assembly of 
triangles in which every triangle shares its sides with each 
of its three neighboring triangles. 

• The area of such a triangulated surface is 

where Si is the area of the i— th triangle and Nt is the total 
number of triangles which compose a given surface. 

• The volume enclosed by this polyhedral surface is the 
summed contribution from Nt tetrahedra 

Nt 

V = Y.y^, 

i = l 

V, = ^S,{njP')i. (4) 

Here Vi is the volume of an individual tetrahedron whose 
base is a triangle on the surface. (rijP-') is the scalar product 
between the outward pointing normal fi to this triangle and 
the mean position vector of the three triangle vertices, for 
which the ji*'' component is given by 

P'^liPl+Pi+Pi). (5) 

The subscript i in Q refers to the i— th tetrahedron, while 
the vectors P\,P2,P3 in 101 define the location of each of 
three triangle vertices defining the base of the tetrahedron 
relative to an (arbitrarily chosen) origin (Fig. Ql. The ver- 
tices are ordered anticlockwise. This ensures that the contri- 
bution to the volume from tetrahedra whose base triangles 

Genus measures the number of handles in excess of the number 
of isolated underdense regions (voids) that the surface of a cluster 
exhibits. 
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lie on opposite sides of the origin add, while the volumes of 
tetrahedra with base triangles lying on the same side of the 
origin subtract out. Thus for both possibilities we get the 
correct value for the enclosed volume (Fig. 0. 

• The extrinsic curvature of a triangulated surface is lo- 
calized in the triangle edges. As a result the integrated mean 
curvature C is determined by the formula 



C 



IE'.- 



(6) 



where £ij is the edge common to adjacent triangles i and j 
and <j)ij is the angle between the normals to these triangles 
hi & hj. 

cos <f)ij = fii.fij. (7) 

The summation in Q is carried out over all pairs of ad- 
jacent triangles. It should be noted that for a completely 
general surface, the extrinsic curvature can be positive at 
some (convex) points and negative at some other (concave) 
points on the surface. To accommodate this fact one asso- 
ciates a number e = ±1 with every triangle edge in Q. e = 1 
if the normals on adjacent triangles diverge away from the 
surface, indicating a locally convex surface, while e = — 1 if 
the normals converge towards each other outside the surface, 
which is indicative of a concave surface. In the former case 
the centre of curvature of the surface lies within the surface- 
body whereas in the latter case the centre of curvature lies 
outside of the surface-body (Fig.|2j. 

• The genus of a triangulated closed polyhedral surface 
is given by the convenient expression 



G = 



X = Nt - Ne + Nv, 



(8) 



where x is the Euler characteristic of the triangulated sur- 
face. Nt , Ne , Nv are, respectively, the total number of tri- 
angles, triangle-edges, and triangle-vertices defining the sur- 
face. 

• A s demonstrated in[S aImi^athvaprakash & Sh andarinI 
il998h : ISathvaprakash. Sahni fc Shandariui 1,19981') ratio's of 
Minkowski functionals called "Shapefinders" provide us with 
an excellent measure of morphology. Therefore, in addition 
to determining MFs we shall also determine the Shapefind- 
ers, T (Thickness), B (Breadth) and L (Length) which are 
defined as follows: 



L = 



C 



(9) 



where G = for simply connected surfaces and G > for 
multiply connected regions. The three Shapefinders (describ- 
ing an object) have dimensions of length and provide us with 
an estimate of the 'extension' of the object along each of 
the three spatial directions: T is the shortest and thus de- 
scribes the characteristic thickness of the object, L is the 
longest and thus characterizes the length of the object; B 
is intermediate and can be associated with the breadth of 
the object. This simple interpretation obviously is relevant 
only for fairly simple shapes. For example, a triaxial ellip- 
soid has the values of T, B and L closely related to the 
lengths of the three principal axes: shortest, intermediate 
and the longest respectively. It is worth stressing that L de- 
fined by eq. @ quantifies the characteristic length between 



holes which should be taken into account when interpreting 
the results (e.g.Tablel^below). 

• An excellent indicator of 'shape' is provided by the di- 
mensionless Shapefinder statistic 



B-T 
B + T' 



L-B 



B' 



(10) 



where P and F are measures of Planarity and Filamen- 
tarity respectively {P,F 5C 1). A sphere has P = _F = 0, 
an ideal filament has P = 0,F = 1 while P = 1,-F = 
for an ideal pancake. Other interesting shapes include 'rib- 
bons' for which P ~ P ~ 1. When combined with the genus 
measure, the triplet {P, F, G} provides an example of shape- 
space which incorporates information about topology as well 
as morphology of superclusters and voids. 

Having presented an ansatz to calculate MFs and 
Shapefinders on a triangulated surface, we now discuss the 
surface generating algorithm (SURFGEN) which creates tri- 
angulated surfaces corresponding to iso-density contours 
evaluated at any desired threshold of the density field. 



3 SURFACE GENERATING ALGORITHM 
(SURFGEN) 

SURFGEN is a versatile and powerful prescription which 
allows us to generate and study surfaces whose physical 
origin can be quite varied and different. Intensity surfaces 
(isophotes), isotherms and iso-density surfaces constructed 
from 3D data provide examples of two dimensional surfaces 
which can be generated and studied using SURFGEN and 
the Minkowski functionals discussed in the previous sec- 
tion*. The present paper is the first of a series in which 
SURFGEN and MFs will be used in conjunction to make 
a detailed morphological study of the supercluster-void net- 
work in the universe. SURFGEN will be applied to a density 
field p(r) defined on a rectangular cubic lattice. 

The SURFGEN algorithm is constructed as follows: 

(i) A point particle distribution (from N-body simula- 
tions, galaxy catalogues) is used to reconstruct the density 
field on the vertices of a cubic lattice via a cloud-in-cell 
(CIC) approach. The lattice itself consists of a large num- 
ber of closely packed elementary cubes. Any given cube is 
characterized by the value of the density at its eight vertices. 

(ii) An appropriate density threshold pth is chosen and 
lattice points at this threshold are found using linear inter- 
polation. The lattice-cells which involve such interpolated 
points along their sides are marked for triangulation which 
is to be done in the next step. 

(iii) SURFGEN determines a closed polyhedral surface at 
the threshold pth by triangulating elementary cubes while 



^ Non-geometrical shape-statistics based on mass moments etc 
can give misleading results whe n applied to large scale structu re 



as demonstrated inlSathvaprak ash. Sahni fc ShandarinI il99d 
* SURFGEN is a modified and extended version of the March- 
ing Cubes Algorithm (MCA) which is used in the field of med- 
ical imaging to render high resolution images of internal or- 
gans by processing the output generated by X-ray tomography 
iLorenson fc Clindll987^ . 
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(A) 



(B) 



Figure 1. Volume enclosed by the surface is estimated by vectorially summing over the volumes of individual tetrahedra; each having 
its base at one of the triangles composing the surface; and its apex at an arbitrarily chosen point which we call the origin. Note that the 
origin can lie both within as well as outside the surface. The left panel shows how SURFGEN estimates the volume when the origin lies 
outside the surface. In this case, the contribution (to the enclosed volume) from triangles falling in the same solid angle carries opposite 
sign for tetrahedra OABC and OFED. Clearly the volume of OFED must be subtracted from the volume OABC in order to give the true 
volume enclosed by the surface. Anticlockwise orientation of the vertices used in determining the normals of the base triangles helps us 
bring this about. The right panel shows the second possibility for which the two triangles ACB and DEF form tetrahedra whose fourth 
(common) vertex lies within the surface. In this case, the contributions from ODEF and OACB add to give the total volume. 




BC = 1 

A A 

cos(|) = Hj. n2 O' 




O 



(a) Convex Neighbourhood (e = 1) 



(b) Concave Neighbourhood (e = -l) 



Figure 2. Illustrated here is the algorithm to compute the local contribution to the Integrated Mean Curvature: (a)AABC and ACBD 
are two triangles sharing a side BC of length £. The normals hi and n2 to the two triangles lie in the planes orthogonal to BC. Their 
projection in the plane of the paper diverges from a point O inside the surface. Thus the triangles correspond to a convex neighborhood 
and e = 1 in Eq. j^J. The angle between the two triangles is denoted as 4>- (b) Shown here are a similar pair of triangles, but with inverted 
sense of the normals, so that the sense of the surface is opposite to that shown in (a). We note that in the present case, the two normals 
converge outside the surface, so that their projections in the plane of the paper would meet at O', which lies outside the surface. Thus 
this represents a concave neighborhood, for which e = — 1 in Eq.J^J. The angle </> represents the internal angle between hi and n2. 



simultaneously moving across the lattice. Triangulated ele- 
mentary cubes are then put together to make up the full 
two-dimensional triangulated surface. 



The process of triangulation is based on the following 
observation. 

For any arbitrary threshold of the density, all elemen- 
tary cubes of the density field fall into one of the following 
three classes: 



(i) those which have all eight vertices below the density 
threshold {underdense cubes - Fig. 

(ii) those which have all eight vertices above the density 
threshold {overdense cubes - Fig.|^II), 

(iii) cubes having both overdense and underdense ver- 
tices {surface cubes - Fig.j^j-III). 

When modelling clusters and superclusters overdense 
cubes will be enclosed by our surface while underdense cubes 
will be excluded by it. (For voids the situation will be re- 
versed.) Contour surfaces at a prespecified density thresh- 
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(I) (II) am 

Figure 3. Given a density field defined on a lattice, the elementary lattice-cubes fall into one of the above three classes: (1) The cubes 
which are completely underdense, (2) the cubes which are completely overdense and finally (3) such cubes which have a few vertices 
overdense and the remaining vertices underdense. Cubes referring to the 3'''* category are triangulated to model the surfaces. 



old will lie entirely within surface cubes. Thus the proper- 
ties of surface cubes are vitally important for this surface- 
reconstruction exercise. 

We now describe how SURFGEN constructs a surface 
at the desired threshold pth by triangulating surface cubes. 
We work under the assumption that the underlying density 
field p{r) is continuous, so that in moving from an overdense 
site to an underdense site (along a cube edge) we invariably 
encounter a point at which p — pTH. SURFGEN classifies 
a given lattice-cube in terms of the number of points (on 
the edges) where p = pth- The precise location of these 
points (found by interpolation) tells us exactly where our 
iso-density surface will intersect the cube. Since there are 
eight vertices and each vertex can either be overdense (1) 
or underdense (0), the number of possible configurations of 
the cube are 2* = 256. Clearly the triangulation must be in- 
variant upon interchanging I's with O's, and this reduces the 
number of independent configurations to 128. Upon invoking 
rotational symmetry this number further reduces to only 14. 
(For instance, although there are eight cubes which have a 
single overdense vertex, any two members of this family are 
related via the three dimensional rotation group. Therefore 
a single scheme of triangulation suffices to describe all mem- 
bers of this group.) These fourteen config urations of triangu- 
lated cubes are shown in Fig. 2] (also see iLorenson fc Clin3 
Jl987l) 1. 

The surface intersecting a given surface cube forms quite 
clearly only a portion of the full iso-density surface which 
we are interested in constructing. By explicitly demanding 
that the density field be continuous across the iso-density 
surface one joins triangles across the faces of neighboring 
cubes, thus constructing the full continuous polyhedral sur- 
face. This then is the basic prescription which may be used 
to construct a closed triangulated surface at any desired 
value of the density threshold. However, a few comments 
about triangulation are in order: (i) the triangulation of the 
cubes having four contrasting vertices is not unique. To illus- 
trate this let us consider Fig.j^where the last cube belonging 
to the middle row of Fig. 2] is reproduced in the left panel. 
We note that this cube is symmetric with respect to an inter- 
change of underdense sites with overdense sites. As a result, 
we could just as well triangulate it according to the prescrip- 
tion shown in the right panel of Fig. |S] This degeneracy in 
the way a given cube can be triangulated is inherent also in 
the last two cubes of the last row of Fig. 2] Noticeably, at 
all such places where we encounter one of these cubes, one 
scheme of triangulation has to be preferred over the other, 
otherwise the surface could become discontinuous. Thus to 




rn (ID 

Figure 5. The cube shown here is the last in the middle row 
of Fig. 4. It has four overdense and underdense vertices and, as 
shown above, it can be triangulated in two distinct ways. Two 
other cubes (not shown) also allow for this possibility (of be- 
ing triangulated in two distinct ways). These are the last two 
cubes belonging to the last row of Fig. 4. The original MCA al- 
gorithm does not mention this degeneracy. However, it can be 
shown that this degeneracy can be broken if we include informa- 
tion about neighboring cubes. Requiring that the density field be 
continuous across the neighboring cubes breaks an occasional de- 
generacy which may arise, and allows one to triangulate surfaces 
unambiguously. 



create a unique triangulated surface one must triangulate 
such cubes tn tandem with their neighbors, ensuring that 
the continuity of the surface is maintained, (ii) There are 
instances when two neighboring cubes have a density con- 
figuration which does not lead to the complete closure of 
the given surface. The resulting surface has a hole (Fig. |S^) 
which must be filled by constructing two additional triangles 
which close the surface (Fig.^t). The probability of hole for- 
mation is finite every time the common face shared by two 
cubes shows two vertices of type 1 along one diagonal and 
of type along the other diagonal. 

Thus, in order to generate a closed polyhedral surface at 
a given value of the density threshold p = pth SURFGEN 
uses the 14 independent triangulations of an elementary lat- 
tice cube in conjunction with an algorithm which ensures 
that the surface is continuous. We should emphasize that 
surfaces generated in this manner need not be simply con- 
nected (see Fig.|7|l. Indeed iso-density surfaces constructed at 
moderate density thresholds in N-body simulations as well 
as 3D galaxy surveys tend to show a large positive value 
for the genus (EqlSJ. This issue will be discussed in detail 
in forthcoming paper(s). To enable an online calculation of 
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Figure 4. Depending upon the value of density (or temperature, if isotherms are being constructed and photon counts, if isophotes are 
being modelled) at its vertices, a surface cube can be triangulated according to one of the fourteen independent configurations shown 
above. The black dots represent overdense vertices whereas underdense vertices are represented by open circles. 





Figure 6. Formation of a 'hole' is illustrated. We 'fill' this hole 
by a pair of triangles as shown in the right panel of the figure. 
Since the original MCA algorithm was meant strictly for visual- 
ization purposes, this problem appears to have been overlooked 
in Lorenson and Cline (1987). We correct for this in our trian- 
gulation scheme by filling holes whenever they occur since its of 
prime importance to work with closed polyhedral surfaces when 
constructing MPs. 



the MPs, our surface modelling scheme also adheres to the 
following requirements: 

• The vertices of all the triangles are stored in an an- 
ticlockwise order. This enforces a uniform prescription on 
normals which always point outward from the surface being 
modelled. Information regarding the directionality of nor- 
mals is of great importance since it is used both for calcu- 
lating the volume as well as the mean curvature. 

• For online calculation of the mean curvature, informa- 



tion concerning a given triangle must be supplemented with 
information concerning the triplet of triangles which are its 
neighbors. This allows us to unambiguously determine the 
contribution to the mean curvature from a local neighbor- 
hood. 

• The total number of triangles, triangle-edges and 
triangle-vertices of a triangulated surface are counted which 
enable us to determine its genus when surface construction 
is complete. 

We might add that, although in this paper the den- 
sity field is reconstructed from N-body particle positions by 
means of a CIC prescription, this is not a necessary prereq- 
uisite for SURFGEN which is versatile enough to be used 
to triangulate iso-density contours of fields reconstructed 
using more elaborate techniques such as Delauney tesse - 
lations (ISchaa p fc Wevgaertll2OO0l : Ivan de Wevgaertll2002D . 
Wiener reconstruction or adaptive smoothing. As it turns 
out, these schemes tend to maintain the structural com- 
plexity in the system without diluting the patterns in an 
isotropic manner. However, since SURFGEN constructs a 
triangu lated surface by 'marchin g across neighboring lattice 
cubes ' llLorenson fc Clind Il987tl the density field must be 
specified on a cubic lattice for this approach to give mean- 
ingful results. ^ 

To summarize, our surface modelling code incorporates 
the triangulation of the entire set of 256 possible configura- 



^ We should emphasize however that the ansatz developed in sec- 
tion |21 for determining the Minkowski functionals for a triangu- 
lated surface can easily be adapted to other triangulation schemes 
such as those described in O'Rourkc 1.1998) ; van dc WcvgaerJ 
i2002l). 
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Figure 7. A cluster appearing in the ACDM simulation is triangulated using SURFGEN. The cluster is multiply connected with two 
tunnels and hence is described by a genus value of two. 

spheres correspond to surfaces of lower constant density in 
this model. We assume io ~ jo ~ ko — 15, for the centre of 
the sphere in our numerical calculations. 

Applying the Shapefinders & HlOjl to a sphere of ra- 
dius R one finds the simple result T = B = L — R and 
F = P = 0. Figure |H] compares the area, volume, mean cur- 
vature and Shapefinders measured for a triangulated sphere 
against known analytical values for these quantities. This 
figure clearly demonstrates that exact values, and values 
computed using triangulation, match exceedingly well down 
to the very lowest scale. 

We should mention that the Genus evaluated using Eq|S| 
is identically zero for spheres of all radii and for all possible 
deformations of an ellipsoid (to be discussed next). This pro- 
vides an excellent independent endorsement of our method- 
ology by demonstrating that our triangulations of a sphere 
and ellipsoid indeed result in closed, continuous surfaces. 

4.2 Triaxial ellipsoid 

A triaxial ellipsoid is an excellent shape with which to test a 
morphological statistic. This is because, depending upon the 
relative scales of the three triaxial ellipsoid can be 

oblate, prolate or spherical. We saw previously that MPs and 
Shapefinders give extremely accurate results for spherical 
surfaces. We now demonstrate that this remains true even 
for surfaces which are highly planar or filamentary. 

The parametric form of an ellipsoid having axes a, b, c 
and volume V = ^ahc is 

r = a(sin 9 cos (j))x + &(sin &sin (f>)y -\- c(cos 9)z (13) 
where 0^0^ 2n, ^ 6 ^ n. For the purposes of our study. 



tions of a cube within a single scheme and makes possible 
online computation of the area, the volume, the integrated 
mean curvature and the genus of a triangulated surface. The 
code has been tested on a variety of standard density dis- 
tributions and eikonal morphologies as well as on Gaussian 
random fields. In the next section, we present results based 
on this analysis. 



4 RESULTS FOR STANDARD EIKONAL 
SURFACES 

In order to test the accuracy of our ansatz we generate trian- 
gulated surfaces whose counterpart continuum surfaces have 
known (analytically calculable) Minkowski functionals. 

4.1 Spherical structures 

In the first exercise, we consider spherically symmetric den- 
sity distributions and generate surfaces of constant density 
for a variety of density-thresholds starting from a chosen 
maximum radius down to grid-size. A convenient density 
distribution for this exercise is 

„rj ^.^ - / 1' ^ (*o,Jo,feo), .... 

P^''^'"' \ po, (i,j,fc) = (io,jo,fco), ^ ^ 

where R is the distance between (i,j,k) and the centre. 

R = V('-^o)2 + (i-jo)2 + (fc-fco)2. (12) 

Since the density field falls ofi^ as the inverse of the dis- 
tance from the centre at (io,jo,fco), any threshold pth is 
associated with a sphere of radius R = Thus larger 
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Figure 8. Minkowski Functionals and dimensional Shapofinders for a spliere. Note that the results for the Minkowski Functionals 
evaluated using our surface modelling scheme (SURFGEN) virtually coincide with the exact values. 



we systematically deform a triaxial ellipsoid and study how 
its shape evolves in the process. When shrunk along a sin- 
gle axis a triaxial ellipsoid becomes planar. Simultaneous 
shrinking along a second axis makes it cigar-like or filamen- 
tary. 



Planarity grows from an initially low value ~ 0.0 to a large 
final value ~ 0.4 as the ellipsoid becomes increasingly oblate. 
Filamentarity, on the other hand remains small at ~ 0.09. 
Both Filamentarity and Planarity are determined to great 
accuracy by SURFGEN. 



4-2.1 Oblate spheroids 

We first study the accuracy of our triangulation scheme 
for planar configurations by considering oblate deforma- 
tions of an ellipsoid. In this case two axes a and b are held 
fixed (a ~ b), while the third axis c is slowly shrunk lead- 
ing to an increasingly planar surface. Our results for this 
case are compiled in two figures. Figure |5] shows MFs as 
they evolve with the dimensionless variable c/a. (The ex- 
act results for the Minkowski functionals which we quote 
ar e based on the analytical express i ons f or MFs given 
in ISahni. Sathvaprakash fc ShandarinI (1 19981) , we refer the 
reader to that paper for more details.) Figure |5| clearly 
demonstrates that the values of MFs obtained using triangu- 
lation match the exact values to a remarkable degree of accu- 
racy. Indeed, for a wide range in c/a the two distinctly differ- 
ent calculational algorithms give virtually indistinguishable 
results, thereby indicating that the triangulation ansatz is, 
for all practical purposes, exact! 

Figure 1101 shows the evolution of all the Shapefinders, 
i.e., T,B,L,P and F together with the percentage errors 
in the estimation of dimensional Shapefinders. We find that 
L is estimated to greatest accuracy with maximum error of 
~ 0.4%, while T and B can be determined to an accuracy 
of ~ 0.8%. We further note (right panels of Figure [TI71 that 



4-2.2 Prolate spheroids 

Next we study prolate deformations of our oblate ellipsoid. 
We start with b a, c <^ a, and shrink the second axis b 
while keeping a and c fixed, so that finally c ~ 6 ^ a and 
our initially oblate ellipsoid becomes prolate. Our results are 
again summarized in two figures. 

Figure 1111 reveals very good agreement between mea- 
sured and true values of MFs, with the former tending to be 
slightly smaller than the latter. 

Turning to the Shapefinders we find that, with the pos- 
sible exception of extremely prolate figures, the Shapefinders 
are remarkably well determined. Indeed, even the extremely 
prolate ellipsoid with axis ratio b/a < 0.2 has a largest error 
in T of only ~ 1.7% while errors in B and L never exceed 
1% (Fig. I12II . This figure also shows the evolution of pla- 
narity (P) and filamentarity {F) as our ellipsoid becomes in- 
creasingly more prolate. In keeping with our expectation F 
steadily increases from a small initial value ~ 0.1 to ~ 0.44. 
Planarity drops from its large initial value to a small final 
value P ~ 0.1 and is slightly underestimated for b/a 0.3. 
We therefore conclude that both P and F are determined 
by SURFGEN to a sufficient accuracy for prolate ellipsoidal 
figures. 
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c/a c/a c/a 

Figure 9. Values of Minkowski functionals determined by the surface modelling scheme SURFGEN are shown along with their exact 
values for oblate deformations of an ellipsoid. Note that for c/a > 0.2, values for Minkowski Functionals evaluated using SURFGEN 
virtually coincide with the exact results. 




Figure 10. Values of Shapcfindcrs derived using the surface modelling scheme SURFGEN are shown along with their exact values for 
oblate deformations of an ellipsoid. The lower middle panel shows percentage errors in the estimation of these quantities. We note that 
the errors are all in the range of ±0.9%. The right panels show the evolution of planarity and filamentarity with continuous shrinking of 
c-axis of the ellipsoid. Note that planarity increases from its low initial value whereas there is only a marginal increase in filamentarity. 



4.3 The Torus and its deformations 

Next we extend our analysis to manifolds which are multiply 
connected by considering the deformations of a torus which 
has an elliptical cross section and which we shall refer to as 
an elliptical torus. A torus is an important surface on which 



to test SURFGEN for two reasons: it is multiply connected 
and, unlike an ellipsoid, it contains regions which are con- 
vex (on its outside) as well as concave (on its inside). The 
elliptical torus can be described by three parameters o, b, c. 
In this case the elliptical toroidal tube has diameter 2nb and 
a & c are its radii of curvature in two mutually orthogonal 
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Figure 11. Values of Minkowski functionals determined by the surface modelling scheme SURFGEN are shown along with their exact 
values for prolate deformations of an ellipsoid. Note that for b/a ^ 0.3, the results for the Minkowski Functionals evaluated using 
SURFGEN virtually coincide with the exact values. 




Figure 12. Values of the dimensional Sliapefinders derived from the surface modelling scheme SURFGEN arc shown along with their 
exact values for prolate deformations of an ellipsoid. The lower middle panel shows percentage errors in the estimation of these quantities. 
We note that the error in the estimation of all the quantities lies within +1% to -2%. The right panels show the evolution of planarity and 
filamentajrity with continous shrinking of the b-axis. Note that the planaxity decreases from its high initial value and the filamentarity 
increases steadily. 



directions. The elliptical torus reduces to the more familiar 
circular torus when a = c. We choose to work with the el- 
liptical torus because changing the values of a, b, c can give 

rise to a large variety of surfaces all of which (by definition) 
are multiply connected but which have very different shapes. 



Thus our surface modelling scheme SURFGEN can be put 
to a rigorous test. The parametric form for the elliptic torus 
is 

r = (6 + csin0) cos9x + (6 + csin0) sinOy + o cos 00, (14) 
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Table 1. Minkowski Functionals for some extreme deformations of an elliptical torus, -^r^, and give the percentage error in 

the determination of the Minkwoski functionals (MFs) using SURFGEN. An accuracy of better than 1% is achieved for the three MFs; 
Curvature (C), Surface area (A) and Volume (V), while the genus is determined exactly. 

b,a,c Surface C ^% S ^% V ^% Genus 

(40.0,37.9,37.9) Sphere-with 7.93 x 10^ +0.47 5.98 x 10'' -0.01 1.13 x lO'' -0.001 1 
-hole 

(140.0,19.9,1.99) Ribbon 2.8 X 10^ +0.03 7.1 x 10* -0.42 1.1 x lO"' -0.830 1 

(60.0,57.9,3.86) Pancake 1.2 X 10^ +0.40 8.8 X 10* -0.18 2.7 x 10^ +0.590 1 

(50.02,3.52,3.52) Filament 9.92 X 10^ +0.52 6.95 x 10^ +0.18 1.21x10* -0.580 1 



Table 2. SURFGEN determined values of the Shapefinders T (Thickness), B (Breadth), L (Length), P (Planarity) and F (Filamentarity) 
describing extreme deformations of an elliptical torus. Also given alongside are the percentage errors in their estimation. 
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-0.48 


31.57 
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0.95 
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Ribbon 


4.67 


1.08 


25.63 


-0.47 


109.96 


0.03 


0.692 


-0.57 


0.622 


0.32 


1 


Pancake 


9.09 


0.74 


74.19 


-0.59 


47.11 


0.41 


0.782 


-0.26 


-0.223 


2.19 


1 


Filament 


5.23 


-0.95 


7.01 


-0.43 


39.49 


0.46 


0.145 


1.40 


0.698 


0.28 


1 



where o, c < 6, ^ < 27r and 6 < 27r. We shall com- 
pare results due to SURFGEN with the exact analytical 
results for four kinds of tori - a nearly spheroidal torus, 
a ribbon, a pancake, and a filament. Tables and |5| re- 
fe r to these four surfaces; relevant figures can be found 
in lSahni. Sathyaprakash fc S handarin flQQSh. TableQshows 
the estimated values of the Minkowski Functionals and the 
percentage error in their estimation. Table |5| shows all the 
Shapefinders for these surfaces^. This table also shows the 
percentage errors in the estimation of the Shapefinders. 

We should emphasize that the genus estimated by SUR- 
FGEN for all these deformations has the correct value of 
unity. This gives an important independent check on the 
self-consistency of the triangulation scheme since, the ab- 
sence of even a single triangle out of several thousand, would 
create an artificial 'hole' in our surface and give the wrong 
value for the genus - a situation that has never occurred for 
any of the several dozen deformations of either the ellipsoid 
or the torus. 

Our results clearly demonstrate that the surface mod- 
elling scheme SURFGEN provides values for the Minkowski 
functionals which are in excellent agreement with exact ana- 
lytical formulae. We have further tested SURFGEN against 
a set of Gaussian random field samples. Owing to a subtle 
issue of the boundary conditions involved in the correspond- 
ing test, we present the corresponding results in an appendix 
at the end of this paper. For the time being, it suffices to 
make a note that our ansatz for studying iso-density con- 



^ It should be noted that the definition of the third 
Shapefinder L is modified fr om that reported in 
ISahni. Sathyaprakash &: ShandarinI il998l) . New definition 
of L is as given in Eq|5] The newly adopted definition of L incor- 
porates treatment of simply connected and multiply connected 
surfaces in within a single scheme. 



tours and their morphology works exceedingly well both for 
simply connected as well as multiply connected surfaces and 
on Gaussian random fields. 

In the next section, we apply SURFGEN to cosmolog- 
ical N-body simulations and give an example of the study 
of the morphology in three models of structure formation: 
ACDM, rCDM and SCDM. 



5 APPLICATIONS: MORPHOLOGY OF THE 
VIRGO SIMULATIONS 

As a first application of the SURFGEN code, we report here 
a morphological study of cosmological simulations. We use 
dark matter distributions of three cosmological models sim- 
ulated by the Virgo consortium. The models are (1) a fiat 
model with Qq ~ 0.3 (ACDM), and two models with fl — 
1, namely (2) one with the standard CDM power spectrum 
(SCDM) and (3) a flat model with the same power spectrum 
as the Qo = 0.3 ACDM model (also referred to as rCDM). 
The shape parameter F (^fioh) for SCDM is 0.5, whereas 
for the other two models F = 0.21. Amplitude of the power 
spectrum in all models is set so as to reproduce the observed 
abundance of rich galaxy clusters at the present epoch. A 
detailed discussion of the cosmological parameters and sim- 
ulations can be found in Jenkins ot al. ( 1998). 

The data consist of 256'^ particles in a box of size 239.5 
h~^Mpc. In order to carry out a detailed morphological 
study, we first need to smooth the data and recover the 
underlying density field. The issue of density field recon- 
struction has been well explored by several groups studying 
the topology of the large scale structure. Here we fol low the 
smoothing technique used by ISpringel et al.l (^9^ which 
they adopted for their preliminary topological analysis of 
the Virgo simulations. We fit a 128^ grid on to the box. 
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Thus, the size of each cell is 1.875 /i~^Mpc. Since there are 
on an average 8 particles per cell, the smoothing is not too 
sensitive to shot-noise. The samples are fairly rich and we 
can smooth the fields in two steps. In the first, we apply a 
Cloud in Cell (CIC) smoothing technique to construct a den- 
sity field on the grid. Next we smooth this field with a Gaus- 
sian kernel which offers us an extra smoothing length-scale. 
There are various criteria to fix this smoothing scale which 
depend on the relative sparseness/richness of the samples. 
To avoid discreteness effect (so that there are no regions in 
the survey where density is undefined) , the simplest criteria 
for sparse samples has been to set A ^ 2.5£g, where £g is the 
grid-size. There are other sophisticated criteria to deal with 
sparse samples which are tested and reported in literature, 
which need not be reiterated here on account of our samples 
being rich. We are left with two choices; either to smooth 
the fields at a scale comparable to ro, the correlation scale, 
or to probe the morphology at smaller length scales with an 
interest to study the regularity in the behaviour of the MPs; 
of which the latter is more promising. The average interpar- 
ticle distance in the present case is ~ 0.94 h~^Mpc, which is 
far outweighed by the grid-size and the correlation length. 
Hence, we adopt a modest scale of smoothing, A = 2 h^^ 
Mpc, which is small compared to usual standards as well as 
enough to provide the necessary smoothing. The Gaussian 
kernel for smoothing that we adopt here is 

In this section we first study the global MFs for all 
the three models. Next we investigate the morphology of 
the percolating supercluster network. This is followed by a 
statistical study of the morphology of smaller structures. We 
summarize our results in the next section. 

5.1 Global Minkowski Functionals 

We scan the density fields at 100 values of the density thresh- 
old pTH, all equispaced in the filling factor 

FFv = j ®iP- P™)d\, (16) 

where Q{x) is the Heaviside Theta function. FFy measures 
the volume-fraction in regions which satisfy the 'cluster' cri- 
terion Pciustcr ^ Pth at a given density threshold pth. In the 
following, we use FFy as a parameter to label the density 
contours. 

^ Since the Gaussian kernel that we use is radially symmetric, 
it could, in principle, diminish the true extent of anisotropy 
in filaments and pancakes. This efi'ect could be minimized by 
considering anisotropic kernels and/or smoothing techniques 
based on the wavelet transform. An even more ambitious ap- 

i iroach is to reconstruct density fields using Delaunay tesselations 
van de Wevgaeri2002^ . Density fields reconstructed in this man- 
ner appear to preserve anisotropic features and may therefore 
have some advantage over conventio nal 'cloud-in-cell' technique s 
followed by an isotropic smoothing <SchaaD fc Weygaerd [20001) . 
Thus, as far as the goal is to utilize the geometry of the patterns to 
discriminate between the models, such smoothing schemes should 
prove more powerful. We hope to return to these issues in forth- 
coming paper(s). 



At each level of the density field (labelled by FFy), 
we construct a catalogue of clusters (overdense regions) 
based on a Friends of Friends (FOF) algorithm.* Next we 
(i) run the SURFGEN code on each of these clusters to 
model surfaces for each of them and (ii) determine the 
Minkowski Functionals (MFs) for all clusters at the given 
threshold (these are referred to in the literature as par- 
tial MFs). Global MFs are partial MFs summed over all 
clusters. Thus, at each level of the density, we first com- 
pute the partial MFs and then the global MFs. Our plots 
of global MFs as functions of FFy are shown in Figures 
1131 This figure shows some interesting features. The three 
cosmological models have appreciably different morphology 
and hence can be distinguished from one another on the 
basis of morph ological measures. Similarly to the perco- 
lation statistic JPominik fc Shandarii]|l992h the difference 
between models is brought out much better if, instead of 
plotting the MFs against the volume fraction (equivalently 
'volume filling factor'), we choose the 'mass fraction' defined 
as FFm = ^^uAai f P Pth) dV . (The 'mass fraction' can 
be thought of as the 'volume or mass filling factor' ^ in (ini- 
tial) Lagrangian space, while the 'volume fraction' is the 
filling factor in (final) Eulerian space.) By employing this 
parameter for studying the global MFs, we essentially probe 
the morphology of the iso-density contours (which may re- 
fer to different thresholds of density or density contrast, but) 
which enclose the same fraction of the total mass. Since clus- 
tering tends to pack up the mass into progressively smaller 
regions of space, such a study connects to aspects of grav- 
itational clustering in a direct way. Thus, we see that at 
all the values of FFm the volume that encloses the same 
fraction of the total mass is least for ACDM and most for 
rCDM. A further advantage of employing this parameter is 
that the behaviour of the rest of the global MFs is no longer 
restricted to follow the same pattern (like it is in case of 
studying global MFs with the volume filling factor as the 
labelling parameter). Thus, not only the peak, but also the 
peak-position and the shape of the MF-curves are sensitive 
to the models being investigated. Thus, the mass fraction 
parametrisation could be more useful in discriminating the 
models from one another as well as comparing the mod- 
els with the observations. Before we commence our analysis 
of individual objects (connected overdense/underdense re- 
gions) we would like to draw the attention of the reader to 
an issue in nomenclature. The 'clusters' and 'superclusters' 
referred to further in this paper are elements constituting 
the 'cosmic-web' and defined on the basis of isodensity con- 
tours in N-body simulations. They should not be mistaken 
for the clusters and superclusters of galaxies seen in the sky. 



* The clusters/superclusters discussed in the present paper are 
defined as connected overdense regions lying above a prescribed 
density threshold. Due to the large smoothing scale adopted, the 
overdensity in clusters ranges from 5 ~ 1 to 5 ~ 10, which makes 
them more extended ( > few Mpc) and less dense than the galaxy 
clusters in (for instance) the Abell catalogue. 

^ Since the density in the Lagrangian space is uniform there is 
no difference between the volume and mass fractions. 
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Figure 13. Global Minkowski functionals for ACDM, rCDM and SCDM are shown as functions of the volume fraction (= volume filling 
factor, FFy). The volume-curve is the same for all models and grows linearly with the volume fraction; this is simply a restatement 
of the definition of 'filling factor'. Notice that the amplitude of the remaining three MFs (Area, Curvature and Genus) is substantially 
greater in SCDM than in ACDM; with rCDM falling between the two. This could mean that large scale structure is much more 'spongy' 
in SCDM with percolating structures in this model showing many more 'holes' (or tunnels) and resulting in a large value for the genus. 
A relative shift towards left, of the position of the peak of the genus curve in ACDM from FFy =0.5 is indicative of the bubble shift, 
implying more amount of clumpyness in this model compared to other two models. It should be noted that for a Gaussian random field, 
the peak occurs at FFy = 0.5. Gravitational clustering modulates the genus curve by lowering its peak and producing a shift to the 
left (the 'bubble shift') or to right (the 'meatball shift') depending upon the other model parameters. The relatively large value of the 
surface area and the curvature may indicate that moderately overdense superclusters have many more 'twists and turns' in SCDM than 
in either tCDM or ACDM, at identical values of the filling factor. In all models the genus curve has a large negative value at large values 
of the volume fraction. The reason for this is that the percolating supercluster, at extremely low density thresholds (high FFy), occupies 
most of the volume. Voids exist as small isolated bubbles in this vast supercluster and lead to a large negative genus whose value is of 
the same order as the total number of voids. In the opposite case, at very high density thresholds associated with small filling fractions, 
we probe the morphology of isolated and simply connected clusters. Since the genus for these clusters vanishes, the global topology of 
large scale structure at high thresholds of the density generically approaches zero for all the models. The vertical lines in the lower right 
panel show the percolation threshold in SCDM (dashed), rCDM (dot-dashed), and ACDM (solid) which occurs at intermediate density 
thresholds corresponding to density contrast <5 1. It is clear from this figure that the percolating supercluster has a relatively simple 
topology at the onset of percolation. The rapid increase in genus value as one moves to lower density thresholds (larger FFy) reflects 
the progressive increase in the 'sponge-like' topology of the percolating supercluster which is more marked in SCDM than in ACDM. 



5.2 Cluster abundance and percolation 

In this subsection we first discuss the percolation properties 
of the three density fields being studied. In this context we 
study how the total number of clusters and the fractional 
volume in the largest cluster vary as we scan through a set 
of density levels corresponding to equispaced fractions of 
volume and total mass. Figure lTsl shows the volume fraction 
and mass fraction as a function of the density contrast. It is 



interesting to note that SCDM and rCDM show the same 
pattern of behaviour which differs from ACDM. It is further 
to note that at high density thresholds {S >, 1) the filling 
factor (at a constant density threshold) is greater in ACDM 
than in SCDM or tCDM, while exactly the reverse is true 
for underdense regions. This figure, which relates FFm,v 
to 5, serves as a reference point to all of our subsequent 
percolation studies. Figure ^1 also shows that, relative to 
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Figure 14. Global Minkowski functionals for ACDM, rCDM and SCDM are shown as functions of the mass fraction defined as 
^^^M = ^^total I P Pth) dV. We see that some of the degeneracy between the different models seen in Fig. 1131 is broken when 
one plots the MFs in terms of the mass fraction -F-Fm instead of the volume fraction FFy. The markers on the three curves show the 
mass fraction at the onset of percolation in each of the three cosmological models. (For the sake of greater clarity, vertical lines through 
the markers are shown in the lower right panel.) 



other models, more mass occupies less space in ACDM. Thus 
almost 67% of the total mass in the ACDM universe resides 
in just 25% of the volume. (In case of rCDM and SCDM 
~ 62% mass occupies ~ 34% volume.) 

A statistical pair which help quantify the geometry of 
large scale structure and its morphology are the 'Number 
of Clusters Statistic' (NCS) and the 'Largest Cluster Statis- 
tic' (LCS); both are shown in Figure \Wl The Number of 
Clusters Statistic shows cluster abundance as a function of 
the density contrast (lower left), mass fraction (lower mid- 
dle) and volume fraction (lower right). The Largest Cluster 
Statistic measures the fractional cluster volume occupied by 
the largest cluster: 

^CS^^, (17) 

where Vi is the volume of the i— th cluster in the sample and 
Vlc is the volume occupied by the largest cluster. Summa- 
tion is over all clusters evaluated at a given threshold of the 
density and includes the largest cluster. In Figure 1161 LCS 
is shown as a function of the density contrast (upper left). 



mass fraction (upper middle) and volume fraction (upper 
right). 

The square, triangle and circle in NCS denote crit- 
ical values of the parameters S = 5ciuster-max, FFy = 
^'^ciustcr-max and F Fm at which the cluster abundance 
peaks in a give model. Similar etchings on the LCS curves 
denote values of 5 = (5pcrc, FFy = FFpac and FFm at which 
the largest cluster first spans across the box in at least one 
direction. This is commonly referred to as the 'percolation 
threshold'. It is important to note that the number of clus- 
ters in SCDM is considerably greater than the number of 
clusters in ACDM at most density thresholds. It is also in- 
teresting that percolation takes place at higher values of the 
density contrast (and correspondingly lower values of the 
volume filling factor) in the case of ACDM. (5porc — 2.3 
for ACDM, 5pcrc ^ 1 for SCDM and rCDM.) Furthermore, 
while in the case of ACDM 5ciustor-max ^ 5perc, in the case 
of SCDM and rCDM, (5ciustcr-max ^ ^perc. Thus in the latter 
two models, as the density threshold is lowered, the cluster 
abundance initially peaks, then, as the threshold is lowered 
further, neighboring clusters merge to form the percolating 
supercluster. In ACDM on the other hand, the threshold at 
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Figure 15. Here we study the evolution of the volume fraction and the mass fraction w.r.t. the density contrast for all the three models. 
ACDM shows maximum density contrast. SCDM and rCDM form another class of models which have similar density contrasts but 
smaller than ACDM. The vertical line across both the panels refers to the mean density, which can be used as a marker to study the 
mass fraction which is contained in the overdense volume (the volume occupied by regions above the mean density). We note that the 
overdense volume in ACDM is '^25% which is about 10% smaller than that occupied by the other two models. At the same time, the 
mass that this volume encloses is ~67% and is about 5% larger than that enclosed by the overdense volume in the other two models. 



which the cluster abundance peaks also signals the formation 
of the percolating supercluster. Having said this we would 
like to add a word of caution: the analysis of one realiza- 
tion alone does not allow us to asses reliably the statistical 
fluctuations in the threshold of percolation. Therefore the 
final conclusion can be drawn only after a study of many 
realizations of each model. 

Clearly both the LCS threshold Sp^ic and the NCS 
threshold (5ciuster-max coutaiu important information and 
most of our subsequent description of supercluster morphol- 
ogy will be carried out at one of these two thresholds. Table|3 
shows the Minkowski functionals and associated Shapefind- 
ers for the 10 largest (most voluminous) superclusters com- 
piled at the percolation threshold for the three cosmological 
models ACDM, rCDM and SCDM while Tablegllists values 
for the Planarity and Filamentarity for these superclusters 
along with the mass that these enclose and their genus value. 

Our (239.5h"^Mpc)^ ACDM universe contains 1334 
clusters & superclusters at the percolation threshold. Of 
these the smallest is quasi-spherical with a radius of a few 
Megaparsec while the largest is extremely filamentary and 



percolates through the entire simulation box (see figure lTTIl . 
From Figures ITTlfclTSl we find that the percolating ACDM 
supercluster is a slim but massive object.^" It contains 4.5% 
of the total mass in the universe yet occupies only 0.6% 
of the total volume. In SCDM the percolating supercluster 
contains 4.4% of the total mass and occupies 1.2% of the 
total volume, i.e., a similar mass occupies two times greater 
volume than in the ACDM. The ~ 10^ ACDM clusters with 
S ^ Spcic occupy 4.4% of the total volume and contain 33% 
of the total mass in the universe. Gravitational clustering 
thus ensures that most of the mass in the ACDM universe 
is distributed in coherent filamentary regions which occupy 
very little volume but contain much of the mass. (Of the 



Surface visualization is a difficult task especially if we wish to 
follow tunnels through superclusters to check whether our visual 
impression of the genus agrees with that calculated using the 
Euler formula jSJ. MATLAB has been used for surface plotting 
and the 'reality' of tunnels is ascertained by rotating surfaces and 
by viewing the supercluster at various angles. 
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Figure 16. Number of Clusters statistic (NCS) and the Largest Cluster Statistic (LCS) are studied as functions of the density contrast 
S (left panels), mass fraction FFm (middle panels) and volume fraction, (equivalently volume filling factor), FFy (right panels). The 
markers on the three curves, in the case of NCS, show values of <5clustcr-max (left), ^'^M.max (middle) and -F-Fv.max (right) at which the 
cluster abundance peaks. In the case of LCS, the markers indicate percolation values of <5pcrc (left), i^i^M.pcrc (middle) and -F-Fv.pcrc 
(right) at which the largest cluster percolates by running through the entire simulation box. 



1334 ACDM clusters and superclusters identified at the per- 
colation tiireshold, the ten most voluminous are extremely 
filamentary and contain close to 40% of the total cluster 
mass (see table 01 . The remaining 60% of the mass is dis- 
tributed amongst other 1324 objects !) 

The percolating supercluster is more filamentary and 
less planar in ACDM [(P, F)ACi3M = (0. 13,0.81)] than in 
SCDM [(P,F)scDM=(0. 14,0.78)]. The percolating super- 
cluster of the rCDM is least filamentary and most planar 
with (P, _F)tCdm=(0.15,0.70). It is interesting that super- 
clusters evaluated at 5porc are more 'porous' (i.e. have a 
larger genus) in SCDM than their counterparts in ACDM. 
For instance the percolating supercluster in ACDM has only 
6 'holes' which tunnel all the way through it, while 20 tun- 
nels pass through the percolating supercluster in SCDM; see 
Table El (This number increases as the density threshold is 
lowered as demonstrated in Figure 1181 ) Another interest- 
ing feature of gravitational clustering appears to be that 
although the volume occupied by clusters at 5ciustcr-max 
does not vary much between models, the fractional mass 
contained in ACDM clusters is almost twice that in either 
rCDM or SCDM; see Figure ITCl 



5.3 Cluster Morphology 

Cosmic density fields contain a wealth of information. As 
demonstrated in Fig. 1161 at density thresholds significantly 
lower than iJporc most clusters merge to form a single perco- 
lating supercluster. On the other hand, the slightly larger 
NCS density contrast, ^dustor-max, provides an excellent 
threshold at which to study the morphology of individual 
objects since it is precisely at (5ciustcr~max that the cluster 
abundance peaks. We study the morphology and topology of 
large scale structure in a two-fold manner: (i) properties of 
all superclusters are analysed at one of the two thresholds: 
5pcrc, 5ciustcr-max; (ii) the morphology of the largest super- 
cluster is extensively probed as a function of the density 
contrast. 

Figure 1181 shows the evolution of the morphology and 
topology of the largest cluster as we scan through different 
threshold levels of the ACDM density field. The density con- 
trast, volume fraction and genus for a few of these levels have 
been labelled and are given in the right corner of the figure. 
We note that at high density thresholds the largest cluster 
occupies a very small fraction of the total volume and its 
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Figure 17. The largest (percolating) supercluster in ACDM. This cluster is selected at the density threshold which marks the onset 
of percolation ((5perc = 2.3). As demonstrated in the figure, the cluster at this threshold percolates through the entire length of the 
simulation box. It is important to note that the percolating supercluster occupies only a small fraction of the total volume and its 
volume fraction (filling factor) is only 0.6%. Our percolating supercluster is a multiply connected and highly filamentary object. Its 
visual appearance is accurately refiected in the value of the Shapefinder diagnostic assigned to this supercluster: (T, B, L) = (5.63, 7.30, 
70.03) h-'^Mpc and (P,F,G) = (0.13, 0.81, 6). 



shape is characterized by significant filamentarity (~0.67) 
and negligible planarity (~ 0). We note a sharp increase 
in the filamentarity of the largest cluster as we approach 
the percolation threshold (shown in the figure as a solid 
triangle). At smaller values of the density contrast, 5 <, 3, 
supercluster filamentarity rapidly drops while its planarity 
considerably increases. The drop in filamentarity of the su- 
percluster is accompanied by a growth in its complexity with 
the result that, at moderately low thresholds {5 ~ 0.5), the 
percolating supercluster can contain several hundred tun- 
nels. At very low thresholds {5 < 0.5) the percolating su- 
percluster is an isotropic object possessing negligible values 
of both planarity (P) and filamentarity (F). 

Since 5porc and (5ciustor-max contain information pertain- 
ing to morphology and connectivity, we compile partial MPs 
for individual clusters at both these thresholds for the three 
cosmological models ACDM, SCDM and rCDM. We find it 
convenient to work with the cumulative probability function 
(CPF) which we define as the normalized count of clusters 
with the value of a quantity Q to be greater than q at a 
given value 'q'. The quantity Q could be one of the MPs or 
one of the Shapefinders. We shall study the dependence of 
CPF with Q on a log-log scale. 

Our results are presented in Figures ITUl and l2Ul for clus- 
ters compiled at the NCS threshold (5ciustGr-max. For large 
values of the Minkowski functionals (V, A, C) the CPP de- 
clines rapidly. Although the curves may look similar, some 
of them are statistically very different from their peers. 
The first four columns of Table |1| show the results of the 



Kolmogorov-Smirnov test applied to the distribution of the 
MPs. (We do not carry out a similar exercise for the genus 
since the vast majority of clusters at the NCS threshold 
are simply connected.) In particular the CPP of the masses 
clearly distinguish the ACDM model from SCDM & rCDM, 
in agreement with Pig. 1191 

Figure I2UI shows the Cumulative Probability Function 
for planarity and filamentarity in our cluster sample. The 
last two columns of Table 0| show the results of the KS test 
for these statistics. For all three pairs they are among three 
best discriminators of the models. The value of KS statistic 
d suggests that the signal comes from relatively small 
values of P (< 0.1) and F (< 0.25). The more conspicuous 
differences in the tail of the distribution is not statistically 
significant due to poor statistics. Clusters in ACDM are the 
least anisotropic, a fact which corresponds to their relatively 
early formation and therefore longer evolution. From Figure 
I20l we also find that clusters at the NCS threshold are signif- 
icantly more filamentary than they are planar. This appears 
to be a generic predic tion of gravitational clustering 

as demonstrated bv [A rnol'd. Sh andarin fc Zel'dovichl 

J1982D: iKlvpin fc Shandarin! ' Jl983l) : 

Se^hyaprakash^ Sahn^^^I^ n^arinl il996l) and 

Bon^Koftnanir^ogos^^'^g^ 

Figure im is a scatter plot of Shapefinders T,B,L for 
clusters in ACDM defined at the percolation threshold. The 
strong correlation between T and B in the left panel in- 
dicates that two (of three) dimensions defining any given 
cluster assume similar values and are of the same order as 
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Table 3. The ten most voluminous superclusters (determined at the percolaton threshold) are listed with the mass they enclose, 
their associated Minkowski functionals (Volume, area, curvature and genus) and Shapefinders T, B, L, P and F. The first row in each 
cosmological model describes the percolating supercluster and appears in boldface. It is interesting that in all three cosmological models 
the top ten superclusters contain roughly 40% of the total mass in overdense regions with 5 (Spcrc • (The precise numbers are 40% for 
ACDM, 37.8% for SCDM and 45.6% in rCDM.) It should also be noted that the mass of a typical supercluster in ACDM is somewhat 
smaller than that in the other two models mainly on account of the fact that the adopted particle mass in ACDM simulations is 
smaller than in simulat ions of SCDM and rCDM. (The particle mass in ACDM is 6.86 x lO^O/i-^iWe and that in SCDM and tCDM 
is 2.27 X 10^^ h~^MQ Ijenkins et al.iri998l) ). It should be noted that the interpretation of L as the 'linear length' of a supercluster 
can be misleading for the case of superclusters having a large genus. In this case L X (G + 1) provides a more realistic estimate of 
supercluster length since it allows for its numerous twists and turns. The morphology of the objects is conveyed through their planarity 
P and filamentarity F. We note that the most voluminous and massive structures are highly filamentary in all the models. The largest 
supercluster is most filamentary in case of ACDM and least so in case of rCDM. 
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the correlation length. (Note that T ~ B ~ 5h~^ Mpc for 
the largest superclusters in Table El) The clustering of ob- 
jects near T ~ _B (P ~ 0) in this panel suggests that our 
clusters/superclusters are either quasi-spherical or filamen- 
tary. The scatter plot between B & L in the right panel of 
Figure 1211 breaks the degeneracy between spheres and fila- 
ments. The mass-dependence of morphology is highlighted 
in this panel in which larger dots denote more massive ob- 
jects. This figure clearly reveals that more massive clus- 
ters/superclusters are, as a rule, also more filamentary, while 
smaller, less massive objects, are more nearly spherical. The 
concentration of points near the 'edge' of the scatter plot 
in the right panel, arises due to an abundance of low mass 



compact quasi-spherical objects with L ~ B (_F ~ 0). The 
large number of massive superclusters with F > 0.5 serves 
to highlight the important fact that the larger and more 
massive elements of the supercluster chain consist of highly 
elongated filaments as much as ~ 100h~* Mpc in length, 
with mean diameter ~ 5h~^ Mpc (see Table |^. It is im- 
portant to reiterate that almost 40% of the total overdense 
mass resides in the ten largest objects listed in table|21 while 
the remaining 60% is distributed among 1324 clusters. 

We also find short filaments to be generally thinner 
than longer ones in agreement with predictions made by 
the adhesion model with regard to the formation of hierar- 
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Shape-Space trajectory of the Largest Supercluster 
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Figure 18. The morphological evolution of the largest supercluster in ACDM is shown as a series of open squares in shape-space {F, P}. 
Each square corresponds to a different value of the density tliresliold which is progressively lowered from a large initial value (5 ~ 6.9; left 
most square) until the mean density level (5 = 0; lower most square). The legend lists the density contrast, the associated volume fraction 
and the genus of the largest supercluster at six monotonically decreasing values of the density contrast (1 — > 6). At the highest {5 ~ 6.9), 
the largest cluster appears to have a large filamentarity and a small planarity. The solid triangle (labelled 3) refers to the percolation 
threshold at which the largest cluster first spans ax;ross the simulation box. After percolation, the filamentarity rapidly decreases from 
a maximum value of 0.81 to 0.0 as the threshold is lowered from 5 ~ 2 to 5 ~ 0. The decline in filamentarity of the largest cluster is 
accompanied by growth in its complexity as revealed by its genus value, e.g. , the fractional volume occupied by the largest cluster at the 
mean density is ~26% and its genus exceeds a thousand. 



Table 4. The Kolmogorov-Smirnov statistic d (first row for each pair of models) and the probability (second row) that the two data sets 
are drawn from the same distribution. The models ACDM, rCDM and SCDM are compared on the basis of three Minkowski functionals: 
volume (V), area (A), integrated mean curvature (C), as well as mass (M) and the Shapefinders: thickness (T), breadth (B), length (L), 
planarity (P). and filaniontarity (F). Boldface lii,y;hlijilits the ttiroo srnaJIcsf pTX>hahi!iti('s for each pair of models. 
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Figure 19. Cumulative Probability Function (CPF) is shown as a function of total mass in superclusters (lower left) and as a function 
of the three Minkowski functionals volume (V), surface area (A) and mean extrinsic curvature (C). Superclusters are defined to be 
connected overdense regions lying above the NCS threshold ^cluster— max- (The cluster abundance peaks at 5 = ^dustor-max (see figure 
1161 which makes this a convenient threshold at which to probe morphology.) The colour type and line typo hold the same meaning across 
all the panels. 



chical filamentary structure durinR Rravitation al clustering 
jKofman. Pogosvan. Shandarin fc Melotdll993) . 

To furtiier probe the morphology of clusters and super- 
clusters we define the notion of shape-space in Figure 1221 
Shape-space is two dimensional with the planarity (of a clus- 
ter) plotted along the x-axis and its filamentarity along the 
y-axis. (One can also incorporate a third dimension show- 
ing the genus.) The first panel in Figure 1^ is a scatter plot 
of P and F for clusters in ACDM. The sizes of dots in the 
middle panel are proportional to cluster mass. We note that 
the most massive structures are also very filamentary. In the 
right panel we try to relate the shape of the structures with 
their topology by scaling the size of the dots with the genus 
value of clusters having a given morphology (P, F) . As shown 
here, clusters which are multiply connected (larger dots are 
indicative of more complicated topologies) are also more fila- 
mentary. Together, the three panels show us that more mas- 
sive superclusters are frequently very filamentary and often 
also topologically quite complex. We also see that a large 
number of less massive superclusters are simply connected 
and prolate. These structures are a few Mpc across along 



their two shorter dimensions and ~ 20 Mpc along the third, 
and therefore have appreciable filamentarity (_F ^ 0.3). 



6 DISCUSSION AND CONCLUSIONS 

This paper presents a new technique for studying the geo- 
metrical and topological properties of large scale structure 
using the Minkowski functionals (MFs) . Given a density field 
reconstructed by appropriately smoothing a point data set 
consisting (say) of a distribution of galaxies, our ansatz: 
(i) constructs closed polyhedral surfaces of constant den- 
sity corresponding to excursion sets of a density field, using 
for this purpose a surface generating triangulation scheme 
(SURFGEN); (ii) SURFGEN evaluates the Minkowski func- 
tionals (volume, surface area, extrinsic curvature and genus) 
thus providing full morphological and topological informa- 
tion of three dimensional iso-density contour surfaces cor- 
responding to a given data set. Evaluated in this manner, 
the Minkowski functionals can be used to study the proper- 
ties of individual objects lying either above (clusters, super- 
clusters) or below (voids) a given density threshold. They 
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Figure 20. Cumulative Probability Punction (CPF) for filamentarity (_F) and the planarity (P) of superclusters selected at the density 
threshold <5clustcr— max- The colour type and line type hold the same meaning for both the panels. 



can also be used to study the morphological properties of 
the full supercluster-void network at (say) the percolation 
threshold. The ratio's of Minkowski functionals (Shapefind- 
ers) are used to probe the shape of iso-density surfaces which 
sample the distribution of large scale structure at different 
thresholds of the density. (The highest density thresholds 
correspond to galaxies and clusters of galaxies, moderate 
thresholds correspond to superclusters while the lowest den- 
sity thresholds characterize voids.) The performance of our 
ansatz has been tested against both simply and multiply 
connected surfaces such as the sphere, triaxial ellipsoid and 
triaxial torus. These three eikonal bodies can be smoothly 
deformed to give surfaces which are spheroidal, pancake-like 
(oblate) and filament-like (prolate) etc. Analytically known 
values of the Minkowski functionals for these surfaces al- 
low us to test both our surface modelling scheme and our 
evaluation of Minkowski functionals from triangulation. Re- 
markably, we find that volume, area and integrated extrinsic 
curvature are determined to an accuracy of better than 1% 
on all length-scales for both simply and multiply connected 
surfaces. The integrated intrinsic curvature (genus) is eval- 
uated exactly. We also find that SURFGEN is remarkably 
accurate at calculating the MFs for Gaussian random fields 
(see the appendix). Having validated the performance of the 



code against eikonal surfaces and Gaussian random fields, we 
apply our method to cosmological simulations of large scale 
structure performed by the Virgo consortium. We study the 
geometry and topology of large scale structure in three cos- 
mological scenarios - ACDM, rCDM and SCDM. All three 
cosmologies are analysed at the present epoch (z=0) using 
global MFs, partial MFs and Shapefinders. (The redshift 
evolution of geometry and topology will be discussed in a 
companion paper.) Our main conclusions are summarized 
below: 

• Using the Minkowski functionals we show that, like 
other diagnostics of clustering, supercluster morphology too 
is sensitive to the underlying cosmological parameter set 
characterizing our universe. Although the three cosmologi- 
cal models considered by us, ACDM, rCDM and SCDM dis- 
play features which are qualitatively similar, the geometrical 
and morphological properties of these models are sufficiently 
distinct to permit differentiation using MFs. We demon- 
strate that studying Minkowski Functionals using mass- 
parameterization significantly enhances their discriminatory 
power. 

• SURFGEN employs a new contour-based method for 
determining topology. This method evaluates the genus for 
an iso-density contour surface directly, by applying the Euler 
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Figure 21. Scatter plot for the pair of Siiapefinders T, B (left panel) and B, L (right panel) defining the morphology of clus- 
tcrs/supcrclusters in the ACDM model. The strong correlation between T and B in the left panel near the line P = indicates that two 
of the three dimensions defining a cluster are equal and of the same order as the correlation length (~ few Mpc). Judging from the left 
panel we find that clusters/superclusters in ACDM are either quasi-spherical or filamentary (since both satisfy T ~ S => P ~ 0). The 
degeneracy between spheres and filaments is lifted by the right panel which is a mass- weighted scatter plot for the Shapefinders B,L. 
Each dot in this panel refers to a cluster and its area is proportional to the fraction of mass in that cluster. The concentration of points 
near the line ¥ = (B = L) reflects the fact that a large number of smaller clusters are quasi-spherical. The more massive structures, 
on the other hand, tend to be fllamentary and the largest and most massive supercluster has F = 0.81. All objects are determined at 
the percolation threshold. 



formula Eq.||HJ to the triangulated surface, and presents a 
significant improveme nt over earlie r grid-based methods of 
determining topology (lMelottill99(f) . 

• A study of percolation and cluster abundance reveals 
several interesting aspects of the gravitational clustering 
process. For all models, on decreasing the density thresh- 
old to progressively smaller values one reaches the criti- 
cal percolation threshold at which the largest supercluster 
runs through (percolates) the simulation box. Percolation is 
reached at moderate values of the density contrast ranging 
from 5pcrc ^ 2.3 for ACDM to (Spore ^ 1.2 for SCDM. The 
abundance of clusters reaches a maximum value at (or very 
near) the percolation threshold and the percolating super- 
cluster occupies a rather small amount of space in all three 
cosmological models. Thus the fraction of total simulation- 
box volume contained in the percolating supercluster is least 
in ACDM (0.6%) and greatest in SCDM (1.2%). When taken 
together, all overdense objects at the percolation thresh- 
old occupy 4.4% of the total volume in the ACDM model. 
For comparison, the volume fraction in overdense regions 
at the percolation threshold is ^ 16% in an idealised, 
continuous Gaussian random field iShandarin fc ZeldovichI 
[1989.) . This fraction can increase upto ~ 30% for Gaus- 
sian fields generated on a grid (lYess fc Shandariiil Il996t 
ISahni. Sathvaprakash fc Shandari nl Il997^ . The fact that 
clusters & supercluster occupy a very small fraction of the 
total volume appears to be a hallmark of the gravitational 



clustering process which succeeds in placing a large amount 
of mass (~ 30% of the total, in the case of ACDM) in 
a small region of space (~ 4%). The low filling fraction 
of the percolating supercluster in ACDM (0.006) strongly 
suggests that this object is either plana r or filamentary 
iSahni. Sathvaprakash fc ShandarirJll997^ and a definitive 
answer to this issue is provided by the Shapefinder statistic. 

• Shapefinders were introduced to quantify the visual im- 
pression one has of the supercluster-void network of be- 
ing a cosmic web of filaments interspersed with large voids 
iSahni. Sathvaprakash fc ShandarirJll99d) . By applying the 
Shapefinder statistics via SURFGEN to realistic N-body 
simulations we have demonstrated that: (i) Most of the 
mass in the universe is contained in large superclusters 
which are also extremely filamentary; the vast abundance 
of smaller clusters and superclusters tends to be prolate 
or quasi-spherical, (ii) Of the three cosmological models 
the percolating supercluster in ACDM is the most filamen- 
tary [F ~ 0.81) and the supercluster in rCDM the least 
{F ~ 0.7). (iii) The percolating supercluster in ACDM is 
topologically a much simpler object than its counterpart in 
rCDM, with the former having only ~ 6 tunnels compared 
to ~ 19 in the latter. Other differences between the models 
are quantified in Figures IT^ — 1221 We also show that among 
various morphological parameters, the Planarity and Fila- 
mentarity of clusters and superclusters is one of the most 
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Figure 22. Shape-Space for ACDM, Left panel (multiplicity): The dots in this panel have area proportional to the number of clusters 
with a given (binned) value of Filamentarity and Planarity. Center panel (mass): The dots have area proportional to the total mass 
contained in clusters having a given (binned) value of Filamentarity and Planarity. Note that more massive superclusters are also more 
filamentary. Right panel (Genus): The dots have area proportional to the total genus value of superclusters with a given (binned) value of 
Filamentarity and Planarity. This figure demonstrates the correlation between the mass of a supercluster, its shape and its genus. More 
massive superclusters are, as a rule, very filamentary and also topologically multiply connected. As in the previous figure, all objects are 
determined at the percolation threshold. 



powerful statistics to discriminate between the models (see 
Table m. 

To summarize, this paper has demonstrated that Minkowski 
functionals and Shapefinders evaluated using SURFGEN 
provide sensitive probes of the geometry, topology and shape 



of large scale structure and help in distinguishing between 
rival cosmological models. Having established the strength 
and versatility of SURFGEN it would clearly be interest- 
ing to apply it to other important issues, related to those 
addressed in this paper, including: (i) geometry and mor- 
phology of underdense regions (voids); (ii) geometry and 
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morphology of strongly overdense regions (clusters); (iii) 
redshift-space distortions of supercluster/void morphology; 

(iv) the time-evolution of the supercluster-void network; 

(v) analysis of superclusters and voids in fully three di- 
mensional surveys such as 2dFGRS and SDSS; (vi) SUR- 
FGEN is also likely to provide useful insights in other phys- 
ical and astrophysical situations in which matter is dis- 
tributed anisotropically such a s in the interstellar m edium 
jLazarian fc Poeosvanl lll99d) : iBoumis. et."all J2002M . We 
hope to return to some of these subjects in the near future. 
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APPENDIX A: GAUSSIAN RANDOM FIELDS 
AND THE ROLE OF BOUNDARY CONDITIONS 

This appendix demonstrates the great accuracy with which 
SURFGEN determines Minkowski functionals for Gaussian 
random fields (hereafter GRFs). 

Before embarking on our discussion, we briefly describe 
the samples that we use and summarise the analytical results 
against which the performance of SURFGEN is to be tested. 

We work with three realizations of a gaussian random 
field with a power law power spectrum (n = — 1) on a i28^ 
grid. Each realization of the field is smoothed with a gaus- 
sian kernel of length A = 2.5 grid-units. All the fields are 
normalised by the standard deviation. This leaves them with 
zero mean and unit variance. The MFs are evaluated at a set 
of equispaced levels of the density fleld which coincide with 
the parameter f on account of a being unity (pth = vu = v). 
V is used to label the levels and is related to the volume fill- 
ing fraction through the following equation: 



FFv{v) = 



~dt. 



(Al) 



The samples are flnally used for ensemble averaging. 

The Minkowski Functionals (MFs) of a GRF are fully 
specified in terms of a length-scale Ac, 



Ac = 



2^g(0) 



(A2) 



Ac can be analytically derived from a knowledge of the power 
spectrum, ft can also be estimated numerically by evaluating 
the variance ^(0) of the field and the variance ^ (0) of any of 
its first spatial derivatives (for more details see Matsubara 
2003). 

For a G RF in three dimensions the four MF s (per unit 
volume) are llTomita. HJI199(1 lMatsubarall2003l) : 
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where 4'(x) is the error function. 
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Let us now describe the boundary conditions which 
we adopt and which incorporate the periodic nature of the 
GRFs generated using a FFT routine. GRFs defined on a 
grid are, by construction, periodic in nature. As a result, two 
overdense excursion sets visibly separated from each other 
inside the box and touching two opposite sides of the box 
would constitute a single cluster. To make this more clear, 
let us turn to Fig. lAll which shows four clusters (labelled 
1, 2, 3 and 4 respectively) in the left panel. The boundaries 
of these clusters — as drawn by an analog of SURFGEN 
in 2-D — are shown in bold, and we note that the sides of 
the box (bold dashed lines) are incorporated when we define 
the clusters. However, due to the periodicity in the density 
field, these four clusters are in fact one single cluster (right 





















Figure Al. A generic 2— D Gaussian random field with four 
clusters touching the boundaries of the box is shown in the left 
panel. SURFGEN-like algorithm in 2— D would treat these as four 
different contours. However, due to the periodic nature of the 
GRF, these four clusters in fact constitute a single bigger cluster 
and the real boundaries of the cluster are as shown by bold lines 
in the right panel of the figure. The dashed solid lines are interior 
to the cluster and are replaced by 2-dimensional planes in 3— D, 
which should not be considered in estimating the area and mean 
curvature. We take care of this cfli'cct while estimating MFs for 
GRFs. 



panel). We note that the portion of the contours due to the 
boundaries of the box occurs in the interior of the actual 
cluster (shown with a bold dashed line). 

It is easy to extrapolate this situation to 3— D. Note 
that the interpolation method employed by SURFGEN will 
correctly estimate the global volume of the excursion set 
even if we worked with the field as it is. In 2— D this is same 
as estimating the area enclosed by the full cluster without 
having to piece the four "clusters" (labelled by 1, 2, 3 and 4) 
together. Our estimation of global genus will also turn out to 
be exact, and as it can be shown, this will hold true even if 
either of these four constituent clusters exhibits a nontrivial 
topology. On the other hand, if we work with the same set 
of contours for estimating global area and mean curvature, 
we will wrongly be adding an excess contribution due to 
the sides of the box. (In 2— D, this translates to adding the 
contribution of the dashed line segments to the perimeter of 
the contour in the right panel of Figure lXTl . The actual area 
and mean curvature contribution comes from the interface 
between overdense and underdense regions (such as the 2— D 
contours drawn in solid bold line in both the panels) as 
opposed to the artificial boundaries of the clusters (shown 
dashed) . 

We note that incorporating periodic boundary condi- 
tions in order to construct a total contour as shown in the 
right panel of Fig. lAll is nontrivial from the numerical point 
of view. When the excursion set touches all the boundaries 
(as, for example, at a very low thresold of density, below the 
percolation threshold), the corresponding algorithm could 
further run into an infinite loop, which makes it impractical 
for implementation. In the light of this, it would be more ma- 
negeable if we restricted ourselves to using contours which 
are closed at the boundaries of the box (like the ones in the 
left panel of Fig. lAlL but devise a way to exclude the excess 
contribution to surface area and mean curvature from the 
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boundaries of the box. A uniform prescription of this type, 
when apphed to all the appropriate clusters at all the thresh- 
olds of density will lead to correct estimation of all the global 
MFs. This then, is the choice of boundary conditions incor- 
porated in our calculation of a GRF. A futher cautionary 
note concerns the familiar W-shape of the genus-curve. Re- 
call that the genus-curve peaks At v — signifying a sponge- 
like topology of the medium at the mean density threshold. 
Note further that G(v = ±^3)= -2G(0)exp (-§) are the 
two negative minima of the genus curve. In fact for \p\ > 
1, the genus is always negative, and approaches zero from 
below when \v\ ^ 1. For v < —1, this negative genus refers 
to isolated underdense regions (bubbles) in the overdense 
excursion set. Similarly, the negative genus for u > 1 should 
be interpreted as being caused by isolated overdense regions 
(meatballs) in the underdense excursion set. To summarize, 
in order to determine the genus for a GRF we need to esti- 
mate the genus of the overdense excursion set for v SCO and 
the genus of the underdense excursion set for z/ > 0, thereby 
making use of the symmetry property of the GRFs. (Note 
however that in our earlier determination of the genus curve 
for N-body simulations in this paper, we have dealt with 
the topology of overdense regions only, because of which the 
genus curve is asymmetric with respect to v = Q and is 
always positive when v > 1.) 

Having elaborated on the nature of boundary conditions 
used in our evaluation of MF's for a GRF, we now present 
our results. 

Figure E2l shows the global MFs averaged over 3 realiza- 
tions of the density field (solid lines with la error bars) along 
with the exact analytical results (shown dotted). There are 
in all 40 levels for every realization. Figure IX2I clear Iv shows 
the remarkable agreement between exact theoretical results 
JASI l - 1IA6I 1 and numerical estimates obtained using SUR- 
FGEN. We therefore conclude that SURFGEN determines 
the MF's of a GRF to great precision. 

The above discussion demonstrated the important role 
played by boundary conditions in reproducing the analyti- 
cal predictions for GRFs. Having tested the performance of 
SURFGEN against GRFs, we are faced with the following 
choice of boundary conditions when dealing with N— body 
simulations and mock/real galaxy catalogues: (i) we could 
either correct for the boundaries when dealing with the 
overdense regions which encounter the faces of the survey- 
volume, or, (ii) we could avoid doing this and treat all mea- 
surements with a tacit understanding that area and mean 
curvature contain an excess contribution which arises be- 
cause of boundary effects. Of the two options, the former 
is computationally more demanding. In addition, its effect 
is appreciable only significantly after the system has perco- 
lated, i.e., when there are a large number of clusters touching 
the boundaries of the box. The corresponding value of the 
filling factor (FFv > FFpcrc) is useful only when one wishes 
to compare two samples through the trends in their global 
MFs. Since the structural elements of the cosmic web (su- 
perclusters, voids) are usually identified near the percolation 
threshold, and are therefore not very sensitive to boundary 
effects, we have chosen option (ii) in our analysis of N-body 
simulations in this paper. Option (ii) is also more suited for 
dealing with real galaxy catalogues which do not satisfy pe- 
riodic boundary conditions. We should also point out that 
since we compare two samples under the same conditions — 



the same volume and dimensions of the box, same resolution, 
etc.— the excess contribution to the area and mean curva- 
ture caused by box boundaries enters in identical fashion for 
both the samples. We therefore prefer option (ii) to (i) and 
all our simulations in Section 5 are analysed without cor- 
recting for contributions from the boundaries. (It may also 
be noted that the contribution from the boundaries becomes 
less important as we deal with larger survey volumes.) 
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Minkowski FuncUonais for GRF(n=-l)(\=S.5 units) 
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Figure A2. Minkowski Functionals for a GRF (P(k)~k~^), smoothed with a gaussian filter with A=2.5 grid-units. Values of the four 
Minkowski fnctionals determined using SURFGEN are shown as solid lines together with the la scatter. Exact analytical results are 
shown as dotted lines. 



